Superinflation, quintessence, and the avoidance of the initial singularity 
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We consider the dynamics of a spatially flat universe dominated by a self-interacting nonmini- 
mally coupled scalar field. The structure of the phase space and complete phase portraits for the 
conformal coupling case are given. It is shown that the non-minimal coupling modifies drastically 
the dynamics of the universe. New cosmological behaviors are identified, including superinflation 
(H > 0), avoidance of big bang singularities through classical birth of the universe from empty 
Minkowski space, and spontaneous entry into and exit from inflation. The relevance of this model 
to the description of quintessence is discussed. 
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I. INTRODUCTION 



The description of the matter content of the cosmos with a single scalar field is appropriate during important 
epochs of the history of the universe 0. In this article, a dynamical system approach to a self-consistent nonsingular 
cosmological history is presented in the framework of the classical Einstein equations with a nonminimally coupled 
scalar field. The complete structure of the phase portrait and of the dynamical behavior is presented for the case 
of a scalar field conformally coupled to the spacetime curvature and with a quartic self-interaction potential. This 
exhaustive analysis is made possible thanks to a reduction of the dynamics to a two-dimensional manifold embedded 
in the original three-dimensional phase space, a general property show earlier by some of the authors Q for a classical 
scalar field in a spatially flat universe with arbitrary self-interaction potentials and arbitrary coupling to the curvature. 
Solutions with special dynamical interest are identified, namely heteroclinic and homoclinic solutions in the reduced 
two-dimensional phase-space, and their relevance to a possible classical birth of the universe from empty space is 
discussed. We recall that heteroclinic trajectories in a phase space correspond to bounded solutions connecting two 
different fixed points. Typically, they play the role of separatrices, determining regions of the phase space with 
qualitative distinct dynamical behaviors. Homoclinic trajectories, on the other hand, correspond to solutions starting 
and ending at the same fixed point. Their relevance to chaotic motions has been intensively discussed in the literature 
0. Despite the fact that the model presented here has no chaotic behavior, its homoclinic solutions will probably 
mark candidate regions of the phase space for chaotic motions if a small perturbation to the equations is introduced. 

Our model consists in a universe filled with a self interacting non- minimally coupled scalar field. A crucial ingredient 
of the physics of scalar fields in curved spaces is their nonminimal coupling to the scalar curvature R of spacetime, 
which is required by first loop corrections || , by specific particle theories [Q , and by scale-invariance arguments at the 
classical level j|| . It is well known that nonminimal coupling dictates the success or failure of inflationary models [[| ; 
more generally, it turns out to strongly affect the cosmic dynamics, which is qualitatively richer than in the minimally 
coupled case. We show, indeed, that nonminimal coupling leads to new dynamical behaviors, such as a regime that we 
propose to call superinflation (H > 0), which cannot be achieved with minimal coupling [Jt), and spontaneous entry 
into and exit from inflation, with or without a cosmological constant. Spontaneous superinflation provides a classical 
alternative to semiclassical birth of the universe from empty Minkowski space |8]|| , which is impossible with minimal 
coupling. 

In the next section, we review our mode l, re cently proposed in Jig], and give some definitions. The aspect of 



the phase portraits are presented in section III. Section IV is devoted to the asymptotic analyses of some especial 



solutions. The last section contains the concluding remarks. 
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II. THE MODEL 



We consider the nonminimally coupled theory described by the action 



S = lj d A x^—g (- - + g^d^dA -2V + £R^ 2 } . ( 1 ) 



where R denotes the scalar curvature, ip is the scalar field, k = 8ttG (G being Newton's constant), and £ is the 
nonminimal coupling constant. A cosmological constant A, if present, is incorporated in the scalar field potential 
V We use the full conserved scalar field stress-energy tensor 

t„„ = d^d v i> - z (VpVv - g^u) a, 2 ) + ic^ 2 - x - gtlv (d a ^d a ^P - 2V{il>)) (2) 

(where G M „ is the Einstein tensor), thereby avoiding the widespread effective coupling k c s — k (l — n^ip 2 ) 1 in the 
Einstein equations G^ v — kT^. We consider here the dynamics of a spatially flat Friedmann- Robertson- Walker 
universe with line element ds 2 — dr 2 — o, 2 (t) (dx 2 + dy 2 + dz 2 ^. This yields the trace equation R = —k (a — 3p), the 
energy constraint 3H 2 = na (which guarantees that the energy density a > 0), and the Klein-Gordon equation. More 
explicitly, 

6 [l - £ (1 - 6f) Kijj 2 } (H + 2H 2 ^ - k (6f - 1) V> 2 - 4kF + 6k ^^ = , (3) 

^ ip 2 + 6^KHipip - 3H 2 (l - <V 2 ) + kV = , (4) 

i) + 3Hip + f Jty + fL = o . (5) 
dip 

The (time-dependent) equation of state of the ip field, rather than being imposed a priori, follows self-consistently 
from the dynamics. In the trace equation (||), the second derivative %[) that appears in the pressure has been replaced 
by its expression given by the Klein-Gordon equation. Clearly, in the set of equations (||)-(|5|), the sub-system (^)-(^) 
is a closed implicit two-dimensional system for ip and H = a/a (note that this dimensional reduction is not possible 
for spatially curved universes Jll|,^l). After solving these implicit equations for ip and H, one has 



^ = -6m±—y/G(H,7P) , (6) 

ZK 



H 



1 



1 + k£(6£- l)ip 2 



3 (2£ - 1) H 2 + 3£(6£ - 1)(4£ - \)nH 2 i? T C(6^ - lWVU + (1 - 2£) kV{iI>) - it&^fr 

dtp 



where 



— - V(yj) + 3£(6£ - l)ffV 



(7) 



(8) 



Due to the energy constraint (^), the trajectories are restricted to a two-dimensional manifold £ in the three- 
dimensional (H,ip,ip\ original phase space, possibly with "holes" (dynamically forbidden regions) corresponding 

to g(H,ip) < (cf. Eq. (||)). £ is composed of two sheets corresponding to the positive or negative sign in Eq. (|^). 
The two sheets smoothly join on the boundary Q = of the dynamically forbidden region. In the present work, we 
restrict to the potential 

vy>) = —4> 2 - -v 4 - — , (9 

k 4 k 

consisting of a mass term, a quartic self-coupling and, possibly, a cosmological constant term. For consistency with 
previous works [^) we use the symbols a = Km 2 /6 (m being the scalar field mass) and lo = — n 2 A/9. Figures 1 and 2 
present some aspects of £ for the potential (||). However, several of our main results do not depend on the details of 
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FIG. 1. Aspect of the two-dimensional manifold E, embedded in the three-dimensional phase space (ip, H, ip), for the potential 
([]). The upper graph represents the "+" sheet, while the lower one depicts the "-" sheet. In fact, they are not disconnected, 
the two sheets join smoothly on the boundary Q=Q of the dynamically forbidden region, corresponding to the showed "holes". 
The lines of constant Q are drawn on the sheets. 




FIG. 2. Zoom, near the origin, of the two-dimensional manifold E, embedded in the three-dimensional phase space (ip, ip, H), 
for the potential (|]). The two sheets join forming two symmetric cones with their apex at the origin. Only the cone corresponding 
to H > is presented here, and the lines of constant H are drawn. 
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III. PHASE SPACE PORTRAITS 



In the following, for simplicity, we project the dynamics of the phase space onto the (H, ip) plane, but the true 
nature of £ should always be kept in mind. The fixed points of the system (|3|)-(|^) include de Sitter solutions with 
constant scalar field 

_ 3(a 2 - Qu) 2 _ 6(a - 6fr) 
" K (Q - 6£a) ' % ~ K (n - 6£a) 1 ' 

(ft =/= 6a£), and the solutions (H, ip) = ^±^/— 3u>/k, O^j . The fixed points ( pX| ) exist also for u> = 0, due to the presence 



of the matter field ip (in this case, the two points (±^/— 3lJJk, OJ collapse into the Minkowski space fixed point, 
which is at the apex of the cones described in Figure 2). Here, we restrict ourselves to the case of conformal coupling, 

£ = 1/6- 

The function 

L^^) = U 2 + ^^-— ip 2 + V(iP) (11) 
2 4k 

is such that dL/dt = —3Hip 2 along the trajectories. For H > 0, L is a Lyapunov function in a region containing the 
origin; the solutions are then confined by closed lines of constant L, implying asymptotic convergence to the fixed 
points on the H axis (from Eq. (Q), if ip and ip vanish, H goes to V(0)). This behavior is confirmed by exhaustive 
numerical simulations |l(J and reported in the following. We first exclude a cosmological constant by setting lu = 0. 
The phase portrait qualitatively differs according to the ratio Q/a. 

A. The case Q = 2a 

The Minkowski space (H, ip, ipj — (0, 0, 0) is a fixed point, attractive for H > and repulsive for H < 0; the projec- 
tions of the de Sitter spaces (±ff , ±^oi 0) are saddle points, i.e. they possess attractive and repulsive eigendirections 
in the phase space (Fig. 3a). They are of two kinds: expanding (Hip > 0) or contracting (Hip < 0). The following 
solutions, present only in this particular case, 

H(t) = y^tanh(^Cr) , tf> = ±ip = ±J^ (12) 

(where C = H + 2H 2 = —R/6 is constant), correspond to heteroclinic straight lines connecting de Sitter fixed points, 
starting along the repulsive eigendirection of one of them and ending along the attractive eigendirection of the other 
(Fig. 3a). They are tangent to the boundary of the forbidden regions at (H, ip) — (0, ±ipo). For \H\ > y/C/2, another 
straight line solution is obtained from the general form 

where w\ and w 2 are integration constants. The nonsingular solutions ( |l2| ) connect a contracting (r — > — oo) de Sitter 
regime to a minimum nonvanishing value of the scale factor (r = 0), and then to an expanding de Sitter regime 
(t — * +oo). 

In addition to these straight lines, we found numerically other heteroclinic solutions: one starting at (H, ip) — (0, 0) 
and ending at (—Ho,ipo), and another one from (Ho,ipo) to (0,0). A third solution starts at the expanding de Sitter 
point and goes to infinity, while another one comes from infinity and arrives to the contracting de Sitter point. The 
phase portrait is symmetric about the origin. 

Near the fixed point (0, 0, 0), numerical analysis confirms the peculiar behavior suggested by the Lyapunov function: 
orbits approaching this point with positive H are attracted to it, bouncing back and forth infinitely many times off 

the Q = boundary in the (H,ip) projection (Fig. 3a). In the space (^H,ip,ipj these orbits are seen to spiral down 

on a cone towards its apex at the origin. The cone results from the union of the two sheets in the vicinity of the 
origin, as it is shown in Figure 3. Along the spiral, the orbit passes almost periodically from one sheet to the other 
with period Tbounce = 27r/m (rbounce is obtained from the asymptotic analysis of the next section. Typically, after a 
few bounces, the period coincides with Tbounce with good accuracy) . A similar behavior for f2 < was reported in the 
earlier numerical analysis of tyj, but using the effective coupling k c s(t) and the variables ip and ip. 



4 




FIG. 3. Qualitative phase portraits for the system (|3j)-(|5|). Shadowed regions correspond to the dynamically forbidden regions 
(G(H,ip) < 0, cf. Eq. (g)). Figures a-e, were obtained by using Q./a, respectively, 2, 5, 3/2, 3/2, and 1/2, and lu = 0. Figure f 
corresponds to the case u> = 1/10 and Q/(a — uo) = 3/2. 

In the H < half-plane, the situation is reversed: orbits starting with H < are repelled by the origin and depart 
from it bouncing off the Q = boundary. 

B. The case Q, > 2a 

The situation (Fig. 3b) is analogous to the previous one, but now the straight hctcroclinic solutions are missing, 
and are replaced by the solution starting at the contracting de Sitter fixed point and escaping to infinity, and by the 
solution coming from infinity and arriving to the expanding de Sitter point. The two-sheeted structure of X implies 
that no actual intersections occur between different orbits in Fig. 3b, which live in different sheets but are projected 
on the same plane. 

C. The case a < Q < 2a 

As shown in Fig. 3c, there are no straight heteroclinic lines but new interesting features emerge. A new heteroclinic 
solution appears starting from the origin and ending in the expanding de Sitter fixed point. As in the previous case, 
the quadrant ip > 0, H < is obtained from the ip > 0, H > one by reflection about the tp-axis and time-reversal 
in Fig. 3c. 

The crucial feature of this case is the appearance of a dense set of homoclinic solutions (Fig. 3d) departing from the 
origin with negative H and returning to it with positive H, going around the forbidden region. Superinflation plays 
a central role along these orbits; only a regime with H > permits the smooth transition from an initial contracting 
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(H < 0) phase to an expanding one (H > 0). This transition occurs at the nonvanishing minimum of the scale 
factor. The behavior of this family of homoclinics, as well as of the other solutions, is universal: they rapidly converge 
in the spiraling region near the origin, irrespective of initial conditions. Since all these homoclinics originate from 
a neighborhood of the Minkowski fixed point due to its own instability with respect to perturbations with H < 0, 
and come back to that point, due to the stability for perturbations with H > 0, this behavior constitute a classical 
alternative to the previously proposed semiclassical birth of the universe from empty space ||||] ■ 



D. The case < Q < a 



The fixed points ( [l0| ) disappear and the only bounded solutions are the homoclinics associated with the origin (see 
Fig. 3e). This situation is therefore the most favorable for the classical spontaneous exit from empty Minkowski space. 



E. The case uj < 



Analogous results hold if a small cosmological constant is present (see Figs. 3f), with the phase portrait being 
classified according to f2/(a — u>), but the fixed point (0, 0, 0) of the w = case splits into two de Sitter fixed points 
with memory of the previous stability properties. Now, the approximate period between two consecutive bounces is 

''"bounce — I ■ V J 



IV. ASYMPTOTIC ANALYSIS 

For the conformally coupled case, the trace and the Klein Gordon equations ([|) and (||), after a rescaling and the 
explicit substitution of the Hamiltonian constraint (Q), read 

H + 2H 2 - ai? + Quo = 
■ij) + 3Hip + 6(a - u))i> - (fi - a)ip 3 = (15) 

We have special interest in the asymptotic behavior of the solutions of (|l^) in the vicinity of the attractive fixed 
points. Let us start with the case ui = 0. In this case, our region of interest is the neighborhood with H > of the 
origin. We search for asymptotic solutions of the form 

t t 

H{t) = iM + irt± + Q {t -3) (16) 
with fi(t), f2(t) ,<?i(£), and g2(t) bounded for large t. Inserting ( [l6|) in the equations ( fllf ) one has 

\ (fi + 6a/i) + ^ (/£ - 2f{ + 3 9l f[ + 6a/ 2 ) - 0(t~ 3 ) 

\g'i + ^ (.92 - .9i + 2gl ~ «A 2 ) - 0(t- 3 ) (17) 



By demanding the exactness of (17) up to t 3 order, one gets immediately that g\ is a constant and fi—A cos(v / 6a£ + 
5), leading to 

/a + 6a/ 2 = AV6a{3 gi - 2) sin(v / fei + 5), (18) 

which has the general solution 

f 2 = B cos(Vfe< + 0) - A ^ 9l _ 2 ^ (co8(V6at + S)VQat- cos(Vfat) sind) . (19) 



G 



In order to guarantee the boundedness of /2, it is necessary to have g\ = 2/3. From the equation corresponding to 
the t~ 2 term in the trace equation, we have 



2 

g' 2 = A 2 acos 2 (Vfat + 5) - -, (20) 



which is solved by 



.92 = ^ ®_ /'cos(v / 6ai + 8) sin(v / 6ai + S) + S) + ( — 
2y6a ^ ' \ 2 



Again, by requiring the boundedness of g 2 for large t, we get the condition A = 2/(3y/a). We have, finally, the 
following asymptotic solution for the ui — case 



, . . 2 cos v 6crf 



i/(0 = | + O(t- 2 ). (22) 



From ( p2| ) we can obtain the characteristic period Tbounce = 27r/v6a. Also, the asymptotic solution (|22j) implies that 
all solutions ending in the origin will behave as matter dominated universes for large t, corresponding to a(t) oc t 2 ! 3 . 

The case with a small cosmological constant can be treated analogously, taking into account that now the relevant 
fixed point is ip = and H — y— 3u>. The search for asymptotic solutions of the form 

t t z 

H(t)^ g0 (t) + g -^l + 9 -^ + O(t- 3 ) (23) 
with fi(t), /2 (t ) , g${t), gi(t), and g2(t) bounded for large t, as before, give rise to the equations 

\ (fx + 3. 9o /( + 6(q - + i (/^ - 2/f + 3Cffi/J + 50/2 - fld/i) + 6(a - w)/ a ) = 0(<- 3 ), 

+ 2 ff0 + 6cj) + 1 (gi + 40031) + ^ (52 - 9i + 4<?o<?i + 2g 2 - af 2 ) = 0(t~ 3 ). (24) 

Our asymptotic equations now are considerably more complicated, but, nevertheless, we can obtain an unam- 
biguous period Tb oullcc - From (pij), one has that g decreases exponentially to the value \J— 3ui for large t, im- 
plying that /1 converge to the form /j = Aexp(— 3y/—3u>t/2) cos(-\/6a + Su/At + 5), from which one can obtain 
Tbounce = 27r/ y/6a + 'Sco/A. The scale factor, in this case, obeys a(t) oc exp y/—3ujt for large t, as in the de Sitter 
universe. 



V. DISCUSSIONS 

The (H,ip) plane is divided into sectors by the straight lines H = ±yja/2ijj, H = zt^/aip H = ±v / 2a - cor_ 
responding, respectively, to H = 0, a = 0, and pressure p = (the if-axis corresponds to p = <r/3). The lines 
H = ±y / a/2 ip mark the transition between inflationary (a > and H < 0) and superinflationary (H > 0) regimes. 
The lines H = ±y/ail> divide regions corresponding to inflation and to decelerated expansion, while H — ±y/2aijj 
divide regions of positive and negative pressures (Fig. 4). The crucial condition for superinflation to occur is that 
the line H = (or parts of it) should belong to the dynamically accessible region Q > of the (H, ip) plane. This 
implies that, for an arbitrary potential V(ip), superinflation corresponds to ipdV/dip < (this result will be discussed 
in detail elsewhere). For our particular potential (^|), this requires ft > 0. The superinflationary behavior occurs only 
once along each homoclinic and brings the solution from the primordial Minkowski neighborhood to the succession 
of eras corresponding to different equations of state (during each bounce) , towards infinite dilution and equation of 
state p = 0. In fact, asymptotic analysis for t — * +00 and for any value of a and f2 (with to = 0) shows that the scale 
factor a(r) exhibits oscillations of concavity corresponding to accelerated and decelerated epochs. These oscillations 
are damped as r — > +00; in this regime, a(r) oc t 2 / 3 and the universe becomes matter dominated. 
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FIG. 4. The plane (H, ip) for the u> — case and the equation of state for ip. The darkest region corresponds to the 
superinflation regime (H > 0). During the bounces of the ip solution in the region H < 0, its equation of state corresponds 
to, respectively, radiation domination (crossing the //-axis), matter domination (p — 0), and re-acceleration (between Q = 
and a = lines. These oscillations are damped as r — > +oo, the universe becomes matter dominated a(r) oc t 2 ^ 3 and tends to 
infinite dilution. 

While it is not claimed here that the evolution of our universe is modeled by an entire orbit of the system (||)-(|5|) 
on the accessible manifold S, the application to specific eras of the cosmological history is intriguing. Indeed, in 
the bounces reported above (during which H < 0), one encounters, respectively, radiation domination crossing the 
ff-axis, matter domination [p = 0), acceleration (a possible quintessence model?) until the next bounce in the (H, ip) 
projection, where this sequence is reversed. 

If we identify one period as our cosmological history, then the reported accelerated expansion of the universe today 
p3[ suggests to locate our epoch in the sector between the line H — \faip and the Q = boundary. The identification 
of the age of the universe (~ 10 17 s) with Tbounce would then yield the scalar field mass m ~ 10 -13 eV, which is 
suggestive of an axion M o r of an ultralight pseudo-Goldstone boson (that the quintessence field should be very light 
was already suggested jQ). 

As a step towards a realistic model, one can include a second scalar field coupled to ip, which has the meaning of 
a fundamental field (as is done, e.g., in hybrid inflation), or mimics a baryonic or other fluid. In spite of the higher 
dimensionality of the phase space, many of the features exposed here for a single scalar field survive [|l5||. Although 
very simple, we think our classical model opens interesting avenues to the understanding of quintessence in terms of 
a non-minimally coupled scalar field. 

Finally, one should keep in mind that our detailed semi-analytical analyses are possible thanks to the reduction 
of the original three-dimensional system to a two-dimensional one on a smooth manifold S, which, in turn, strongly 
indicates the absence of chaos (confirmed by an exhaustive numerical analysis). Generically, such a reduction is 
not possible for the case of a spatially curved spacetime and/or several matter fields non-minimally coupled to the 
curvature, despite of many similar results. This is object of present investigations pa]. 
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